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Abstract 

Very recent inelastic a-scattering data on the isoscalar monopole and dipole strength distribu- 
tions in ^^Fe, ^^Ni, and ^°Ni are analyzed in the relativistic quasiparticle random-phase approx- 
imation (RQRPA) with the DD-ME2 effective nuclear interaction (nuclear matter compression 
modulus Kjim = 251 MeV). In all three nuclei the calculation nicely reproduces the observed 
asymmetric shapes of the monopole strength, and the bimodal structure of the dipole strength 
distributions. The calculated centroid and mean energies are in very good qualitative agreement 
with the experimental values both for the monopole, and for the low- and high-energy components 
of the dipole transition strengths. It is noted, however, that while DD-ME2 reproduces in detail 
the excitation energies of the giant monopole resonances (GMR) in nuclei with A > 90, the the- 
oretical centroids are systematically above the experimental values in lighter nuclei with A < 60. 
The latter can be reproduced with an effective interaction with a lower value of ^ 230 MeV 
but, because of the asymmetric shapes and pronounced fragmentation of the monopole strength 
distributions, isoscalar GMR data in light nuclei cannot provide accurate estimates of the nuclear 
matter compression modulus. 

PACS numbers: 21.60.Ev, 21.60.Jz, 21.65.-hf, 24.30.Cz 
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Experimental excitation energies of compressional (monopole and dipole) vibrational 
modes in atomic nuclei can in principle be used to deduce the value of the nuclear mat- 
ter compression modulus K„m [l]- This quantity is related to the curvature of the nuclear 
matter equation of state at the saturation point, and controls basic properties of atomic nu- 
clei, the structure of neutron stars, the dynamics of heavy-ion collisions and of supernovae 
explosions. The nuclear matter compressibility cannot be measured directly, but rather de- 
duced from a comparison of experimental excitation energies of isoscalar giant monopole 
resonances (ISGMR), with the corresponding values predicted by microscopic nuclear effec- 
tive interactions characterized by different values of Knm- 

Inelastic a-scattering experiments have been used in high precision studies of the sys- 
tematics of ISGMR in nuclei with A > 90. Nuclear structure models provide a consistent 
description of the shapes of strength distributions and the mass dependence of excitation 
energies, and thus relate the ISGMR to the nuclear compressibility and to the nuclear matter 
compression modulus. There is much less experimental information, and only few micro- 
scopic theoretical analyses of the structure of compressional modes in lighter nuclei with 
A < 90. While in heavy nuclei the shape of the ISGMR strength distribution is typically 
symmetric, for A < 90 the ISGMR display asymmetric shapes with a slower slope on the 
high energy side of the peak, and with a further decrease of the mass number the ISGMR 
strength distributions become strongly fragmented. An interesting question, of course, is 
whether studies of compressional vibrations in lighter nuclei can provide additional infor- 
mation on the nuclear matter compression modulus. Namely, K.„m corresponds to bulk 
nuclear compressibility, whereas one expects that surface compressibility plays an increas- 
ingly important role in the structure of ISGMR in lighter systems. In a very recent study 
P] isoscalar giant resonances in ^^Fe, ^^Ni, and ^°Ni have been studied with small-angle 
inelastic a-scattering. In particular, most of the expected isoscalar EO has been identified 
below 40 MeV excitation energy, and between 56% and 72% of the isoscalar El strength has 
been located in these nuclei. It was noted, however, that there are no specific microscopic 
calculations of EO and El strength distributions in ^^Fe and ^''Ni. The mass dependence 
of the ISGMR excitation energies between A = 40 and A = 90 was thus compared with 
results of leptodermous expansions based on Hartree-Fock + RPA calculations with Skyrme 
interactions Js^, and constrained relativistic mean- field calculations j4|. The purpose of this 
work is to perform fully self-consistent relativistic quasiparticle random-phase approxima- 
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tion (RQRPA) calculations of isoscalar EO and El strength distributions in ^^Fe, ^^Ni, and 
^°Ni, using a modern effective density-dependent interaction which is known to reproduce 
the systematics of compressional modes in heavier nuclei with A > 90. 

Theoretical studies of nuclear compressional modes in the last decade have employed the 
fluid dynamics approach j^, the Hartree-Fock plus random phase approximation (RPA) ^, 3, 

, linear response within a stochastic one- 



body transport theory 111, the relativistic transport approach [12| , and the self-consistent 
relativistic RPA isl H^Jlsl IgI ]. As has been pointed out by Shlomo et al., however, most 
current implementations of the non-relativistic RPA are not self-consistent, and based on 
numerous approximations [17]. Very recent studies have emphasized the importance of 
a fully self-consistent description of ISGMR, and confirmed that the low value of K^^ = 
210—220 MeV, previously obtained with Skyrme functionals, is an artefact of the inconsistent 
implementation of effective interactions j9|, |l8[- The excitation energies of the ISGMR in 
heavy nuclei are thus best described with Skyrme and Gogny effective interactions with 
Knm ~ 235 MeV. In Ref. it has been shown that it is also possible to construct Skyrme 
forces that fit nuclear ground state properties and reproduce ISGMR energies, but with 
Knm ~ 255 MeV. In Ref. ^9,] a new set of Skyrme forces was constructed that spans a 
wider range of values of Knm and the symmetry energy at saturation density 04. RPA 
calculations with these forces have shown that the ISGMR data are best reproduced with 
Knm = 230 — 240 MeV, whereas higher values of Knm would require unrealistically large 
value of 04. On the other hand, it appears that in the relativistic framework the interval 
of allowed values for Knm is more restricted. A recent relativistic RPA analysis based on 
modern effective Lagrangians with explicit density dependence of the meson-nucleon vertex 
functions, has shown that only effective interactions with Knm = 250 — 270 MeV reproduce 
the experimental excitation energies of ISGMR in medium-heavy and heavy nuclei, and that 
Knm ~ 250 MeV represents the lower limit for the nuclear matter compression modulus of 
relativistic mean- field interactions [16]. 

Data on the compressional isoscalar giant dipole resonance (ISGDR) could also be used 



201 ]. The problem, however, is that 



to constrain the range of allowed values of Knm jisl 
the isoscalar El strength distributions display a characteristic bimodal structure with two 
broad components: one in the low-energy region close to the isovector giant dipole resonance 
(IVGDR) 2huj), and the other at higher energy close to the electric octupole resonance 
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3huj). Theoretical analyses have shown that only the high-energy component represents 



compressional vibrations 



22| |. whereas the broad structure in the low-energy 



corresponds to vortical nuclear flow associated with the toroidal dipole moment 



re gion 
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25l |. However, as has also been pointed out in the recent study of the interplay between 



compressional and vortical nuclear currents 



24| , a strong mixing between compressional and 



vorticity vibrations in the isoscalar El states can be expected up to the highest excitation 
energies in the region ^ 3huj. Nevertheless, models which use effective interactions with 
Knm adjusted to ISGMR excitation ener gies in heavy nuclei, also reproduce the structure of 



the high-energy portion of ISGDR data 121, l26L l27|. 



In this work the isoscalar EO and El strength distributions for the open-shell nuclei 
56pg^ ^^Ni, and ^°Ni are calculated in the relativistic quasiparticle random-phase approxima- 
tion (RQRPA), formulated in the canonical single- nucleon basis of the relativistic Hartree- 
Bogoliubov (RHB) model Q]- The RQRPA model is fully self-consistent: the same in- 
teractions, in the particle-hole and particle-particle channels, are used both in the RHB 
equations that determine the canonical quasiparticle basis, and in the RQRPA equations. 
In both channels the same strength parameters of the interactions are used in the RHB 
and RQRPA calculations. This is an essential feature of the RHB+RQRPA approach and 
it ensures that RQRPA amplitudes do not contain spurious components associated with 
the mixing of the nucleon number in the RHB ground state, or with the center-of-mass 
translational motion. 

In Fig.^we display the isoscalar monopole strength distributions for ^^Fe, ^^Ni, and ^''Ni. 
The RHB+RQRPA calculation has been performed with the DD-ME2 effective interaction 
2^ in the particle- hole channel and, as in most applications of the RHB model js^, the 
finite-range Gogny force has been used in the particle-particle channel. DD-ME2 belongs 
to a new class of relativistic effective nuclear interactions with density-dependent meson- 
nucleon vertex functions. In a number of recent studies it has been shown that this type 
of effective interactions provides a realistic description of asymmetric nuclear matter, neu- 
tron matter and finite spherical and deformed nuclei. These interactions allow for a softer 
equation of state of nuclear matter (i.e. lower incompressibility) and a lower value of the 
symmetry energy at saturation. In addition to nuclear matter and ground-state properties 
of spherical nuclei, the parameters of DD-ME2 have been adjusted to the excitation energies 
of the ISGMR and IVGDR in ^°^Pb. For DD-ME2 the nuclear matter compression modulus 
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amounts Knm = 251 MeV. The strength distributions in Fig. 2]can be compared with the 
data from Ref. (Figs. 8, 9 and 10). In all three nuclei the calculation predicts asymmetric 
shapes for the isoscalar EO strength distributions, in agreement with data. In particular, 
an additional tail in the transition strength is obtained above the main ISGMR peaks for 
^^Fe and ^°Ni. For ^^Ni most of the strength is distributed over two major peaks, with an 
additional pronounced high-energy tail. The arrows denote the positions of the experimental 



centroid (^i = mi/mo) and mean energies (i^a = yms/mi), where rrik = f E^R[E)dE are 
the energy moments, and R{E) is the transition strength distribution function. We note 
that in all three nuclei the main ISGMR peak predicted by the RQRPA calculation is located 
in the narrow energy window between the Ei and -^3 experimental energies. 

In the upper panel of Fig. |21 we plot the RHB+RQRPA results for the ISGMR centroid 
energies {mi /mo) of a series of spherical nuclei from ^''Ca to ^°^Pb, calculated with the 
DD-ME2 effective interaction, in comparison with data from the Texas A&M University 
, bi ba, y and Osaka ^ blj compilations. We note that the latter data cor- 
respond to peak energies and, especially in nuclei in which a high-energy tail is found above 
the main peak, these values should be somewhat below the TAMU centroid energies. The 
agreement between the excitation energies calculated with DD-ME2 and the TAMU data 
is remarkable for nuclei with A > 90, whereas the theoretical centroids are systematically 
above the experimental values in lighter nuclei. The origin of this discrepancy is not un- 
derstood, but it could be due to the fact that in light nuclei the surface incompressibility 
plays a more important role in determining the ISGMR, whereas Knm represents the volume 
incompressibility. The former quantity is seldom taken into account when adjusting the pa- 
rameters of an effective interaction and, therefore, we do not really expect that DD-ME2 can 
reproduce in detail the moments of asymmetric and even fragmented isoscalar EO strength 
distributions in light nuclei with A < 60. 

It seems that data on ISGMR in light nuclei are not very useful in extracting information 
on the nuclear matter compression modulus Knm- Nevertheless, we have tried to repro- 
duce these data with few additional effective interactions. In the recent analysis of nuclear 
matter incompressibility in the relativistic mean- field framework families of density- 
dependent interactions with different values of the nuclear matter compression modulus 
Knm and symmetry energy at saturation (volume asymmetry) a4, were adjusted to repro- 
duce nuclear matter and ground-state properties of spherical nuclei. By performing fully 
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consistent RRPA/RQRPA calculations of isoscalar EO and isovector El strength distribu- 
tions in spherical nuclei with A > 90, it has been shown that the comparison with data 
restricts the values of K^^ to ~ 250 — 270 MeV, and the range of volume asymmetry to 
32 MeV < a4 < 36 MeV. A weak correlation between a4 and Knm was found, i.e. inter- 
actions with lower volume asymmetry allow for slightly lower values of Knm- Therefore in 
addition to DD-ME2, the family of interactions with = 32 MeV and K„m=230, 250, and 
270 MeV Q has been used in a RHB+RQRPA calculation of ISGMR in ^°Ca, ^^Fe, ^^Ni, 



^°Ni, and ^°Zr. The resulting mean energies £"3 = \/jn^mi are plotted in the lower panel 
of Fig. 121 in comparison with data from Refs. P, E3, We notice that while DD-ME2 



and the K„m=250 MeV effective interaction reproduce the experimental value for ^°Zr, 
data in lighter nuclei are better described by the effective interaction with K„m=230 MeV, 
except possibly for ^^Ni, but for this nucleus the experimental E3 differs considerably from 
the values in the neighboring ^^Fe and ^°Ni j^. 

For the DD-ME2 effective interaction, the RHB+RQRPA isoscalar dipole transition 
strength distributions in ^^Fe, ^^Ni, and ^°Ni are shown in Fig. El In all three nuclei the 
El strength is strongly fragmented and distrubuted over a wide range of excitation energy 
between 10 MeV and 40 MeV, in agreement with the experimental results of Ref. Q. In 
the experiment between 56% and 72% of the isoscalar El strength has been located in these 
nuclei below 40 MeV excitation energy, and some missing strength probably lies at higher 
energies. Similarly to the results obtained for heavier nuclei 22 , [2^ , the El strength 
is basically concentrated in two broad structures: one in the region 10 MeV < E^,. < 20 
MeV, and the high-energy component above 25 MeV and extending above 40 MeV excita- 
tion energy. Only the high-energy portion of the calculated El strength is sensitive to the 
the nuclear matter cornpression modulus of the effective interaction. In a number of recent 



theoretical studies l23 



2J, it has been shown that the low-lying El strength mostly cor- 



responds to vortical flow (dipole toroidal mode), although a strong mixture of compressional 
and vortical velocity fields is predicted in the intermediate and high-energy region. 

In Fig. ini the thick arrows denote the locations of the experimental centroid energies 
(mi /mo) in the low- and high-energy regions of the isoscalar El strength in ^^Fe, ^^Ni, and 
^°Ni j2i]- These are compared in Fig. ^ with the theoretical values of the centroids of the low- 
and high-energy components, for different values of Ec, the somewhat arbitrary parameter 
which separates the low- and high-energy regions. We notice a good qualitative agreement 
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between the calculated and experimental centroids in the high-energy region, especially 
taking into account that the El strength above = 40 MeV has not been observed in 
the experiment. In the low-energy region, however, the theoretical centroid energies are 
systematically below the experimental values by MeV, depending on the choice of 



Ec- This effect is in agreement with previous RRPA calculations in heavier nuclei [25|, and 
supports the picture of pronounced mixing between compressional and vorticity vibrations 
in the intermediate region of excitation energies. 

In conclusion, we have performed RHB+RQRPA calculations of the isoscalar monopole 
and dipole strength distributions in A?«60 nuclei. The results obtained with the DD-ME2 
effective interaction (nuclear matter compression modulus K„m = 251 MeV) have been com- 
)ared with very recent data on the EO and El strength distribution in ^^Fe, ^^Ni, and ®°Ni 
^. For the isoscalar monopole resonance we find very good qualitative agreement between 
theory and experiment, both for the asymmetric shapes of the strength distributions, as 



well as for centroid {Ei = mi /mo) and mean energies {E3 = yms/mi). It has been noted, 
however, that while there is an excellent agreement between the ISGMR excitation energies 
calculated with DD-ME2 and the data for nuclei with A >90, the theoretical centroids are 
systematically above the experimental values in lighter nuclei with A < 60. Even though 
because of asymmetric shapes and pronounced fragmentation, ISGMR data in light nuclei 
are probably not very useful for extracting information on the nuclear matter compression 
modulus, we have shown that the ISGMR centroids in nuclei with A < 60 are better de- 
scribed with an effective interaction similar to DD-ME2, but with a lower value of Knm ~ 
230 MeV. The isoscalar El strength distributions calculated with DD-ME2 are in good 
agreement with the experimental results P], and reproduce the observed bimodal structure 
with two broad components in the 2huj and 3huj energy regions. The calculated centroid 
energies of the low- and high-energy El components in ^^Fe, ^^Ni, and ^°Ni qualitatively 
reproduce the experimental values obtained from small-angle inelastic a-scattering data. 
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FIG. 1: The RHB+RQRPA isoscalar monopole strength distributions in ''"Fe, ''^Ni, and ''"Ni 
calculated with the density-dependent effective interaction DD-ME2. The experimental centroid 
(mi/mo) and mean {^/m^/mi) energies obtained from (a, a') scattering Q| are denoted by grey 
and black arrows, respectively. 
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FIG. 2: In the mass region 40 < ^ < 208 the RQRPA results for the ISGMR centroid energies 
(m-i/mo), calculated with the relativistic DD-ME2 effective interaction, are plotted as a function of 
mass number and compared with data from the Texas A&M University (TAMU) 
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Q,|3],|32, 



33] and 



27l | compilations (upper panel). In the lower panel the calculated ISGMR excitation 



ener^s i/ms/mi for several medium-mass nuclei are shown in comparison with the data from 
Ref. 2l. In addition to the DD-ME2 interaction, three additional effective interactions with the 

n 

values of the nuclear matter compressibility K„m=230, 250, and 270 MeV [Ig] have been used in 
the RHB+RQRPA calculation of the ISGMR strength distributions. 
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FIG. 3: The RHB+RQRPA isoscalar dipole transition strength in ^''Fe, ^^Ni, and ''"Ni calculated 
with DD-ME2 effective interaction. The arrows denote the positions of the experimental centroid 
energies of the low- and high-energy components Q]. 
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FIG. 4: Comparison of the RHB+RQRPA (circles and diamonds) and experimental (thick lines) |j| 
centroid energies of the low- and high-energy components of the isoscalar dipole transition strength 
in ^^Fe, ^^Ni, and ^'^Ni. The results are plotted for different values of Ec, the parameter which 
separates the low- and high-energy regions. 
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